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Описана детерминированная модель распределения дефектов. Проведён сравнительный анализ теоретиче- 
ского расчёта и экспериментальных данных многократного отражения ультразвуковых волн на системе трёх 
цилиндрических дефектов, находящихся в упругом теле конечных размеров. 
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Введение. При сканировании скопления дефектов в упругих телах раздельно-совмещенными 
ультразвуковыми датчиками [1] информация 06 их местоположении и форме может быть получе- 
на на основе характеристик однократно и многократно переотраженных волн. Полное исследова- 
ние такой задачи возможно только на основе решения трехмерной задачи высокочастотной ди- 
фракции ультразвуковых волн. С учетом применяемых частот в ультразвуковом неразрушающем 
контроле (УЗНК) материалов теоретическое исследование трехмерных задач многократного рас- 
сеяния упругих волн может быть проведено в рамках геометрической теории дифракции [2, 3]. 
Особый интерес представляет исследование задач высокочастотной дифракции в рамках стохас- 
тической модели распределения дефектов в теле конечных размеров. Однако в качестве началь- 
ных вариантов изучения этой проблемы естественным является ее изучение в рамках детермини- 
рованных моделей распределения дефектов канонической формы. 

Детерминированная модель и экспериментальные измерения. Модель реализована в ла- 
бораторном макете, который в качестве тела конечных размеров использует стальной образец в 
форме прямого параллелепипеда, толщина образца 40 мм, высота 78 мм, длина 300 мм. На рас- 
стоянии Н от верхней грани, на которой располагается датчик продольных ультразвуковых волн 
(УЗ) диаметром 30 мм и частотой 2,5 МГц, по дуге окружности радиусом РК просверлены три оди- 
наковых соосных цилиндрических отверстия радиусом г с образующими перпендикулярными 60- 
ковыми стенками. Схема образца с расположением отверстий и датчика приведена на рис.1. Для 
проведения эксперимента используется дефектоскоп-приставка «ЭВУД-ПК» (в дальнейшем де- 
фектоскоп) общего назначения по ГОСТ 23049-94. Он предназначен для УЗ контроля продукции 
на наличие дефектов типа нарушения сплошности и однородности металлов. Дефектоскоп рабо- 
тает совместно с 1ВМ компьютером. 

Лучевое решение задачи. Расположение цилиндрических отверстий в образце и ультразвуко- 
вого датчика продольной волны позволяют измерять амплитуду перемещений в обратно отра- 
женных волнах от поверхностей полостей в продольных волнах, в двукратно отраженных от де- 
фектов 1 и 3 волнах, в трехкратно отраженных волнах от всех трех препятствий. Траектории ука- 
занных лучей приведены на рис.1. Заметим, что эти траектории представляют собой ломаные ли- 
нии, лежащие в одной плоскости. Общие случаи однократного и многократного отражения про- 
дольной волны от поверхностей, находящихся в упругой среде, получены в работах [3, 4]. Явные 
выражения перемещений в отраженных волнах получены в рамках геометрической теории ди- 
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фракции [5] на основе модификации физической теории дифракции Кирхгофа [6] асимптотиче- 


ской оценкой дифракционных интегралов методом многомерной стационарной фазы [7]. 
Пьезоэзлемент 





Рис.1. Схема образца с датчиком и осциллограмма с экрана дефектоскопа: 1 — зондирующий импульс; 
2 — импульс, отражённый от отверстия 1 (импульс, отражённый от отверстия 3 по времени совпадает с импульсом, 
отражённым от отверстия 1); 3 — импульс, отражённый от отверстия 2; 4 — импульс, отражённый от отверстия 1 
и переотражённый отверстием 3; 5 — импульс, отражённый от отверстия 1 и переотражённый отверстиями 2 и 3 


Выпишем эти формулы для продольных перемещений и, в случаях однократного, двукрат- 
ного, трехкратного отражения волны. 
Случай однократного отражения: 


и = ВУ (у) езр и, (4+1 )+ (5:” у :) 








ВР 
= (1-Е) +2 (+1, ) (22, 05? у +? у) соз ' у+4Р К; (1) 
2 — № — У у РЁ = У Хх | , 
В=О,К (4лик?) 


Здесь Хо - точка, в которой помещен источник продольной волны, вызванной единичной силой 


—14 Пе 
О.е ‚ Х -точка приема волны; Г. — с, ых, — с, С „>С, - вОоЛНОВЫе числа И сКО- 
рости соответственно продольной и поперечных волн; (@) - частота колебаний; 2 - модуль 
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сдвига; И (У )- коэффициент /?— р отражения (продольной волны в продольную); 


К =А К, - гауссова кривизна; М, = ( К +К, )/ 2 — средняя кривизна поверхности в 


—/ 


*ж 
точке зеркального отражения У : А - кривизна нормального сечения поверхности отра- 
* 


жателя плоскостью луча Хо —У —Х . Кривизна К определяется формулой Эйлера 


о р 25 _ 05а с. „д св 
К = К, с0$ ф +, ЗШ °ф (605 ф=— эй Ф=— —__) выражающей кривизну 
и 7 т 7 


произвольного нормального сечения через главные кривизны А1, А> и угол Ф , образуе- 
мый касательной этого нормального сечения с первым главным направлением. 


Формула (1) дает главный член асимптотики продольных перемещений при К, рт >> Г, 
К, вЫ К, ВЕ, К, К, >> [, где Ки К. - главные радиусы кривизны поверхности в 
точке зеркального отражения у. 

Угол 7 - это угол между направлением 4’ — {—с05@,-с0$ 2,-с0$ у} падения и от- 


хх 
ражения продольной волны и нормалью к поверхности в точке у в локальной декартовой сис- 
теме координат. 





— 


2 & 
Рис.2. М-кратное отражение продольной волны вдоль луча № — У У —... — Ух - Хи 
от граничных поверхностей Х/ полостей, находящихся в упругой среде 


Случай двукратного (№=2) и трехкратного (№=3) отражения (рис.2): 


м 
Л (хр 
О Им ехр| /! р, +1 (5% +2м) 
ис ‚) = Иры: НЙ и | _ С0$ | и= 
р 4лиК | )с037, П | 


и=0 


(2) 





‚4е (ру) 
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Здесь о = — ет - разность между числом положительных и отрицательных собственных 


ь В»... ВЕНЕ: 
значений матрицы Гессе у = (4, 1 =1.2.....2№ (М№=2, М=3), которая является сим- 


метричной со следующими ненулевыми элементами 4,1 < ] (приводятся для случая 
М№=3): 
диагональные элементы: 


ди (е кче.) 
и (ви...) 286.) 


И. > 
внедиагональные элементы: 


—0 - —0 — 
Е уе г у 
Я. 12. — бы - (2 и— те 2) (4-7, (4. А . 


|4. ы —0 ри = 
а з я = - 
Я =. от = “ Ч, 1 . мт} АЙ 
|4. - [ч, - =. и | 2 
Я. 1.242 -- (4, > Г (4 2Уиы | о (7 пт } : а. 
—0 — —0 р — 
Зы Чи» Л а Чи 2 И = (1 > У и | ; И = | й / 


ао —- И в > и. вы 2 а — ре 2 ы } а п- 1, _ Е 


Здесь А”, А” (и=1,2,3) — главные кривизны поверхностей, а /», /„,К» — орты локальной декар- 
товой системы координат, определяемой касательными к главным линиям кривизны и нор- 

% 
малям к поверхностям дефектов в точках зеркального отражения у, и = 1, 2,3 (см. рис.2). 


—0 


Координаты ортов заданы в некоторой глобальной декартовой системе координат, 4, - на- 
ва 
правление падения волны в точке у,. 


Выше выписаны элементы матрицы Гессе И: ’ в случае трехкратного отражения волны. 


Для случая двукратного отражения элементы матрицы Гессе В.” — (4, (1,)=1,2,3,4) находятся 


в левом верхнем углу матрицы />?. 
Получим из общих формул (1) и (2) явные выражения для амплитуды перемещений, при- 
нимаемой датчиком (см. рис.1) продольной волны. Во всех случаях 
И — р = 1) _ ОО) а: > 
К О 
1. Обратное однократное отражение (см. рис.1): 
У=0,Н, =0,5,/„=- 
а) от каждой из полостей 1 и 3 (первый отраженный импульс): 


в г 


А о м о. (3) 
0 ] 2 о В я 





6) от полости 2 (второй отраженный импульс): 
венки = Вин | 
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2. Двукратное отражение волны от поверхностей двух цилиндров 1 и 3 в плоскости, пер- 
пендикулярной их осям (третий отраженный импульс, рис.1): 
д К 2 3 
и =ЯУ, = —, =. =Н+—-Г 2_1 ь 23. К+ 58 2 | : 
4 | - 2 © й я, я. 


и@з) И =В-0,5/.01)/| 6 де | (4) 


3. Трехкратное отражение волны от поверхностей трех цилиндров 1, 2 и 3 в плоскости, 
перпендикулярной их осям (четвертый отраженный импульс, рис.1): 
Иа К 
м. =ЁЬ=В; Е 


= В.0, 5/07) (Н/К+0,5) Ю\/3 | с(с+1(с+5)(<-7) || 


Сравнительный анализ. Измерения проведены на стальном образце. Скорости продольных и 
поперечных волн в стали соответственно равны с, =5850 м/с, с, =3260 м/с. При частоте, излу- 


7 — У> == 


чаемой продольной волной 2,5 МГц, ее длина составляет 2,34 мм. Приведем, например, сравне- 
ние теоретических расчетов и результатов экспериментальных измерений амплитуды в обратно 
отраженных волнах от полостей 1 и 3 и в двукратно отраженной волне последовательно от де- 
фектов 1 и 3. 

В результате эксперимента получены следующие значения амплитуды перемещений в от- 


раженных волнах: амплитуда перемещений и — в обратной однократно отраженной волне от де- 


ю Я ю 
фекта Г (см. рис.1) составляет 60,5 ДБ; амплитуда перемещении и. ’в двукратно отраженной 
волне последовательно от дефектов 1 и П составляет 29,5 дБ. Разность между этими значениями 

| 1.3 
составляет 31 дБ. При этих данных отношение и‘ а. РЗ, 


МА 


Теоретические расчеты амплитуд и’ и и - основывались на аналитических выражени- 


С) 
ях (3) и (4). При этом для отношения и, Дис получено выражение 


снесен 


ее. а = 
[н+%- РЯ `) 

), 
1, =48,94 мм, Г. =12,88 мм, Н +0,5(К +") =50 мм, Г, =0,36 


Е Ре р=2/2 Г. /#. 

Расхождение для отношений амплитуд перемещений в обратно отраженной и двукратно 
отраженной волне экспериментальных измерений и теоретических расчетов составляет 3 дБ. 

С учетом того, что стальной образец изготовлен с погрешностями и излучаемая датчиком 

волна не является плоской, а отношение длины волны (2,34 мм) к диаметру полости (3 мм) не 
является малым, что не в полной мере соответствует лучевой теории, полученный результат яв- 
ляется вполне удовлетворительным. 
Заключение. Используемые в УЗНК основные частоты зондирующих импульсов 2,5-10 МГц по- 
зволяют применять для тел ограниченных размеров аналитические выражения геометрической 
теории дифракции высокочастотных упругих волн, полученные для бесконечной среды, содержа- 
щей препятствия. 
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